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In juxtaposition with the standard model of rotation powered pulsar, the model of vibration 
powered magnetar undergoing quake-induced torsional Alfven vibrations in its own ultra strong 
magnetic field experiencing decay is considered. The presented line of argument suggests that 
gradual decrease of frequencies (lengthening of periods) of long-periodic pulsed radiation detected 
from a set of X-ray sources can be attributed to magnetic-field-decay induced energy conversion 
from seismic vibrations to magneto-dipole radiation of quaking magnetar. 



I. INTRODUCTION 

There is a common recognition today that the standard 
(lighthouse) model of inclined rotator, lying at the base 
of current understanding of radio pulsars, faces serious 
difficulties in explaining the long-periodic (2 < P < 12 
s) pulsed radiation of soft gamma repeaters (SGRs) and 
anamalous X-ray pulsars (AXPs). The persistent X-ray 
luminosity of AXP/SGR sources (10 34 < L x < 10 36 erg 
s _1 ) is appreciably (10-100 times) larger than expected 
from a neutron star deriving radiation power from the en- 
ergy of rotation with frequency of detected pulses. Such 
an understanding has come soon after the detection on 
March 5, 1979 of the first 0.2 second long gamma burst [l| 
which followed by a 200-seconds emission that showed a 
clear 8-second pulsation period [2| and association of this 
event to a supernova remnant known as N49 in the Large 
Magellanic Cloud Q. This object is very young (only a 
few thousand years old) but period of pulsating emission 
is typical of a much older neutron star. In works [3, Q it 
has been proposed that discovered object, today desig- 
nated SGR 0526-66, is a vibrating neutron star, that is, 
the detected for the first time long-periodic pulses owe 
their origin to the neutron star vibrations, rather than 
rotation as is the case radio pulsars. During the follow- 
ing decades, the study of these objects has been guided 
by idea [1, 0] that electromagnetic activity of magnetars 
both AXP's and SGR's, is primarily determined by decay 
of ultra strong magnetic field (10 14 < B < 10 16 G) and 
that a highly intensive gamma bursts are manifestation 
of magnetar quakes [a-flOl]. 

In this work we investigate in some details the model of 
vibration powered magnetar which is in line with the cur- 
rent treatment of quasi-periodic oscillations of outburst 
luminosity of soft gamma repeaters as being produced 
by Lorentz-force-driven torsional seismic vibrations trig- 
gered by quake. As an extension of this point of view, in 
this paper we focus on impact of the magnetic field de- 
cay on Alfven vibrations and magneto-dipole radiation 
generated by such vibrations. Before so doing, it seems 
appropriate to recall a seminal paper of Woltjer 11| who 



was first to observe that magnetic-flux-conserving core- 
collapse supernova can produce a neutron star with the 
above magnetic field intensity of typical magnetar. Based 
on this observation, Hoyle, Narlikar and Wheeler [l2[ 
proposed that a strongly magnetized neutron star can 
generate magneto-dipole radiation powered by energy of 
hydromagnetic, Alfven, vibrations stored in the star af- 
ter its birth in supernova event (see, also, [l3|). Some 
peculiarities of this mechanism of vibration powered ra- 
diation has been scrutinized in our recent work [l4j |. de- 
voted to radiative activity of pulsating magnetic white 
dwarfs, in which it was found that necessary condition for 
the energy conversion from Alfven vibrations into elec- 
tromagnetic radiation is the decay of magnetic field. As 
was stressed, the magnetic field decay is one of the most 
conspicuous features distinguishing magnetars from nor- 
mal rotation powered pulsars. It seems not implausible, 
therefore, to expect that at least some of currently mon- 
itoring AXP/SGR - like sources are magnatars deriving 
power of pulsating magneto-dipole radiation from energy 
of Alfvenic vibrations of highly conducting matter in the 
ultra strong magnetic field experiencing decay. 

In approaching Alfven vibrations of neutron star in its 
own time-evolving magnetic field, we rely on the results of 
recent investigations |15l - [l8| of both even parity poloidal 
and odd parity toroidal (according to Chandrasekhar ter- 
minology (HI) node- free Alfven vibrations of magnetars 
in constant-in-time magnetic field. The extensive re- 
view of earlier investigations of standing-wave regime of 
Alfvenic stellar vibrations can be found in [20j . The spec- 
tral formula for discrete frequencies of both poloidal and 
toroidal a- modes in a neutron star with mass M, radius 
R and magnetic field of typical magnetar, B14 — B/10 14 
G, reads 
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where numerical factor sg is unique to each specific shape 
of magnetic field frozen in the neutron star of one and 
the same mass M and radius R. 
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II. ALFVEN VIBRATIONS OF MAGNETAR IN 
TIME-EVOLVING MAGNETIC FIELD 



III. VIBRATION POWERED RADIATION OF 
QUAKING MAGNETAR 



In above cited work it was shown that Lorentz-force- 
driven shear node-free vibrations of magnetar in its 
own magnetic field field B can be properly described in 
terms of material displacements u obeying equation of 
magneto-solid-mechanics 

pii(r,i) = -L[V x [V x [u(r,f) x B(r,t)}]] x B(r,tj(3) 
u(r,t) = [w(r,t)xr], w(r, t) = A t [V X {r)] d(t). (4) 

The field u(r,t) is identical to that for torsion node-free 
vibrations restored by Hooke's force of elastic stresses 
[H [H| with x(r) = A e f e (r) Pi (cos 6) where f t (r) is the 
nodeless function of distance from the star center and 
Pi(cos9) is Legendre polynomial of degree I specifying 
the overtone of toroidal mode. In the amplitude a(t) 
is the basic dynamical variable describing time evolu- 
tion of vibrations which is different for each specific over- 
tone; in what follows we confine our analysis to solely one 
quadrupole overtone. The central to the subject of our 
study is the following representation of the time-evolving 
internal magnetic field 



B(r,t) = B(t)b(r), 



(5) 



where B{t) is the time-dependent intensity and b(r) is di- 
mensionless vector-function of the field distribution over 
the star volume. Scalar product of (J} with the separable 
form of material displacements 



u(r, t) = a(r) a(t) 



(6) 



followed by integration over the star volume leads to 
equation for amplitude a(t) having the form of equation 
of oscillator with time-depended spring constant 



Ma{t) + K{t)a(t) = 0. 



(7) 



The total vibration energy and frequency are given by 
Ma 2 {t) K,{B{t))a 2 {t) 



E A {t) 



W(t) - \ ^ = B(t)K, K 
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It follows 



dE A (t) a 2 (t) dK(B) dB(t) 



dt 



dB 



dt 



(8) 
(9) 

(10) 



This shows that the variation in time of magnetic field 
intensity in quaking magnetar causes variation in the vi- 
bration energy. In next section we focus on conversion of 
energy of Lorentz-force-driven seismic vibrations of mag- 
netar into the energy of magneto-dipole radiation. 



The point of departure in the study of vibration-energy 
powered magneto-dipole emission of the star (whose ra- 
diation power, V, is given by Larmor's formula) is the 
equation 

=-*>(*), V (t) = ^{t). (11) 

We consider a model of quaking neutron star whose 
torsional magneto-mechanical oscillations are accompa- 
nied by fluctuations of total magnetic moment preserv- 
ing its initial (in seismically quiescent state) direction: 
fi = /i n = constant. The total magnetic dipole moment 
should execute oscillations with frequency equal to 
that for magneto-mechanical vibrations of stellar matter 
which are described by equation for a(t). This means 
that Sfj,(t) and a(t) must obey equations of similar form, 
namely 

(12) 
(13) 

It is easy to see that equations (|T2"j) and ([Tl?)) can be 
reconciled if 



6ji(t) +uj 2 (t)5fi{t) = 0, 

&(t) + u 2 (t)a(t) = 0, co 2 (t) = B 2 (t)n 2 



SfJb{t) = fia(t). 



(14) 



Given this, we arrive at the following law of magnetic 
field decay 

2n 2 n 2 , . 

7 = — - — ^ = constant, (15) 



dB{t) 



dt 
B{t) = 



= -lB 3 (t), 



3Mc 3 
= 2 7 S 2 (0). 



(16) 



The last equation shows the lifetime of quake-induced vi- 
brations in question substantially depends upon the in- 
tensity of initial (before quake) magnetic field B(0): the 
larger B(0) the shorter r. For neutron stars with one 
and the same mass M — 1.4M Q and radius R = 15 km, 
and magnetic field of typical pulsar B(0) = 10 12 G, we 
obtain r ~ 3 x 10 7 years, whereas for magnetar with 
B{0) = 10 14 G, t ~ 7 x 10 3 years. 
The equation for vibration amplitude 

a(t) + u 2 (t)a{t) = 0, 
u 2 (0) 



u/(t) = 



l+t/T 



w(0) = CJ^K. 



(17) 
(18) 



with help of substitution s = l + t/T is transformed to 
the equation 

sa"{s) + f3 2 a(s) = 0, (3 2 = uj 2 {Q)t 2 = const. (19) 

permitting general solution [22j. The solution of this 
equation, obeying two conditions a(t — 0) = ao and 
a(t = t) = 0, can be represented in the form 

a(t) =Cs 1 / 2 {J l {z(t)) - v Yi(z(t))}, z(t) = 2f3s 1 / 2 (t) 
JMr)) 



C = a [J 1 {z(0))-7]Y 1 (z(0)} 



(20) 
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FIG. 1. (Color online) The figure illustrates the effect of mag- 
netic field decay on the vibration frequency and amplitude of 
quadrupole toroidal a-mode presented as functions of x — t/r. 



where J\(z) and Y\(z) are Bessel functions [23| and 



= 2E A (0) 

Mw 2 (0) K(0) 



2 2E A (0) 
a - 



u; 2 (0) 



K(0) 
M 



(21) 



Here by Ea{0) is understood the average energy stored 
in torsional Alfven vibrations of magnetar. If all the 
detected energy -Eburst of X-ray outburst goes in the 
quake- induced vibrations, -Bburst = Ea, then the initial 
amplitude cxq is determined unambiguously. The impact 
of magnetic field decay on frequency and amplitude of 
torsional Alfven vibrations in quadrupole overtone is il- 
lustrated in Fig.l, where we plot a(t) with pointed out 
parameters f3 and rj. The magnetic-field-decay induced 
lengthening of period of pulsating radiation (equal to pe- 
riod of vibrations) is described by 



P(t)=P(0)[l + (t/r)}^, P{t) 



P(0) 



2r [l + (i/r)]V2 



P 2 (0) 
2P(t)P(t)' 



P(0) = 



2vr 



,B(0) 



(22) 



On comparing r given by equations (|16|) and (|22l) . one 
finds that interrelation between equilibrium equilibrium 
value of the total magnetic moment [i of a neutron star 
of mass M = 1.4 Mq and radius R — 10 km vibrating in 
quadrupole overtone of toroidal a-mode is given by 



fi = A(M,R)JP(t)P(t), A 



3 Ale 3 



8tt 2 



fx = 3.8 x 10 37 dP(t) P(t), G cm 3 



(23) 
(24) 



For a sake of comparison, in the considered model of 
vibration powered radiation, the equation of magnetic 
field evolution is obtained in similar fashion as that for 



the angular velocity fl(t) does in the standard model of 
rotation powered neutron star which rests on equations 



dE R 
dt 
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F J/4 2 (f), 
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E R (t) = ^IQ 2 (t), I=-MR 2 
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6fi(t) = [n(t) x [n(t) x fj,]], 

fi = fin, fi 
which lead to 



1 i 
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(25) 

(26) 
(27) 
(28) 



K = 



3/c 3 ' 
1 = 2Kfl 2 (0) 



fi± = (i sin9, (29) 
(30) 



where 9 is angle of inclination of fi to fi(t). The 
time evolution of P(t), P(i) and expression for r are 
too described by equations (|2"2"|) . It is these equa- 
tions which lead to widely used exact analytic esti- 
mate of magnetic field on the neutron star pole: B = 

[3/c 3 /(2tt 2 R e )] l l 2 Jp{t) P(t). For a neutron star of 
mass M = Mq, and radius R = 13 km, one has 



B = 3.2 10 19 yP(t) P{t) G. Thus, the substantial phys- 
ical difference between vibration and rotation powered 
neutron star models is that in the former model the elon- 
gation of pulse period is attributed to magnetic field de- 
cay, whereas in the latter one the period lengthening is 
ascribed to the slow down of rotation [24], HE| . 



IV. SUMMARY 

The last two decades have seen a growing understand- 
ing that magnetic field decay is the key property distin- 
guishing magnetars from pulsars. The magnetic fields 
frozen-in the immobile matter of pulsars (most of which 
are fairly stable to starquakes) operate like a passive 
(unaltered in time) promoter of their persistent radia- 
tion along the axis of dipole magnetic moment of the 
star inclined to the rotation axis. The presented treat- 
ment of spin-down effect emphasizes kinematic nature of 
time evolution of pulsar magnetic moment as brought 
about by rapid rotation with gradually decreasing fre- 
quency and showing that rotation does not affect in- 
tensity of its internal magnetic field. Unlike pulsars, 
magnetars (SGR-like sources) are isolated quaking neu- 
tron stars whose seismic and radiative activity is thought 
to be dominated by decay of magnetic field. Perhaps 
the most striking manifest of seismic vibrations of mag- 
netar are quasi-periodic oscillations of outburst lumi- 
nosity rapidly decreasing (in several hundred seconds) 
from about Lx ~ 10 44 erg s _1 in giant flare to about 
Lx ~ 10 erg s _1 in quiescent regime of long-periodic 
X-ray emission. As was emphasized, this lowest value 
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of luminosity is well above of the rate of energy of ro- 
tation with period equal to that of detected emission. 
With all above in mind, we have set up a model of vibra- 
tion powered radiation of quaking magnetar in which the 
key role is attributed to magnetic field decay. The most 
striking outcome of presented line of argument is grad- 
ual decrease of frequencies (lengthening of periods) of 
magneto-dipole radiation (pulsating with the frequency 
of toroidal a-mode of seismic vibrations) owing its origin 
to magnetic field decay. Remarkably that this prediction 
of the model is consistent with data on gradual decrease 
of frequency of pulsed X- ray emission of such magne- 



tars as IE 2259.1+586 [H and XTE J1810-197 [27|. It 
is hoped, therefore, that presented here theoretical re- 
sults can be efficiently utilized as a guide in observational 
quest for vibration powered neutron stars. 
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